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ABSTRACT 

In the standard model of gas giant planet formation, a large solid core (~ 10 times the Earth's mass) 
forms first, then accretes its massive envelope (100 or more Earth masses) of gas. However, inward 
planet migration due to gravitational interaction with the proto-stellar gas disk poses a difficulty 
in this model. Core-sized bodies undergo rapid "Type I" migration; for typical parameters their 
migration timcscale is much shorter than their accretion timescale. How, then, do growing cores avoid 
spiraling into the central star before they ever get the chance to become gas giants? Here, we present 
a simple model of core formation in a gas disk which is viscously evolving. As the disk dissipates, 
accretion and migration timescales eventually become comparable. If this happens while there is still 
enough gas left in the disk to supply a jovian atmosphere, then a window of opportunity for gas giant 
formation opens. We examine under what circumstances this happens, and thus, what predictions 
our model makes about the link between proto-stellar disk properties and the likelihood of forming 
giant planets. 

Subject headings: protoplanetary disks, origin: solar system — planets: formation, planets:extrasolar 
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The observed correlation of planet occurrence and host star metallicity (Gonzalez 1997; Valenti & Fischer 2005), 
£f) ■ together with the lack of evidence (Quillen 2002; Dotter & Chaboyer 2003; Fischer & Valenti 2005) for accretional 
' pollution of stellar convective zones (Laughlin & Adams 1997; Murray et al. 2001) implies that solids in the disk aid 
VO giant planet forrmation. This in turn provides support for the core-accretion model of giant planet formation, wherein 
' a large protoplanet, ~ 10 Earth masses (M®), forms first, then accretes a massive gas envelope (Mizuno et al. 1978; 
Pollack et al. 1996). The growth of such a large body in itself provides a challenge to our understanding of planet 
accretion. It may require some or all of the following: large disk masses (Lissauer 1987; Thommes et al. 2003), local 
enhancements in solids surface density (Stevenson & Luninc 1988; Cuzzi & Zahnlc 2004) and a substantial fraction of 
planetesimals at <C km size in order to increase the effectiveness of aerodynamic gas drag, thus increasing the capture 
+3 . cross-section of the core atmosphere (Inaba et al. 2003) and making planetesimal random velocities very small (Rafikov 
2004). However, a further problem plagues the core accretion model. Gravitational interaction with the parent gas disk 



causes the inward migration of embedded planetary bodies (Ward 1997; Papaloizou & Larwood 2000; Tanaka et al. 
2002); this is referred to as Type I migration. The migration rate increases with planet mass until the planet becomes 
large enough, of order 10 2 M®, to open a gap in the disk (Goldreich & Trcmaine 1980; Lin & Papaloizou 1986), thus 



switching to what is called Type II migration. The migration rate is therefore peaked in the mass range of cores, of 
^ ■ order 10 1 M®, and for most disk models far outstrips the accretion rate. In order to prevent bodies from falling into 
the star long before they reach core mass, studies of concurrent growth and migration of giant planets have thus far 
resorted to assuming a large reduction in the Type I rate (Alibcrt et al. 2005), or setting it to zero (Trilling et al. 2002; 
Ida & Lin 2004). Though arguments can be made for doing so (see §7), it is nevertheless troubling if this is really the 
only way to salvage core accretion. Thus, our purpose here is to address the question of whether unmitigated Type I 
migration really does doom giant planet cores. To do this, we construct a simple semi-analytic model of concurrent 
core accretion and Type I migration in a viscously evolving protostellar disk. We find that although cores cannot 
form in a young disk, the situation changes as the disk dissipates. We conclude that, depending on disk properties, 
there exists a limited timespan, after the disk becomes sufficiently tenuous but before it can no longer supply a jovian 
atmosphere, within which giant planet cores can successfully grow. 

2. GROWTH AND MIGRATION RATES 

Since we want to concurrently calculate the migration and accretional growth of protoplanets, we require analytic 
expressions for both. For the former, we use the result of Tanaka et al. (2002): 

W = (2.7+ 1.10 W ^{-) Or (1) 

for a mass M body orbiting a mass M* star at radius r, embedded in a gas disk with surface density S s and sound 
speed c s . fl is the Kcplerian angular velocity and (3 = — rflogE 9 /c?logr. For the latter, we assume "oligarchic" growth 
(Kokubo & Ida 1998, 2000) and make use of the corresponding growth rate estimate of Thommes et al. (2003). For 
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simplicity, we neglect the radial motion of planctcsimals due to aerodynamic gas drag, the enhanced capture radius of 
growing protoplanets due to formation of an atmosphere, and fragmentation of planetesimals. The growth rate is then 



DM 



Dt 

where 



AT, m M 2 / 3 , (2) 



A _ 50 (C D b Pg )V^M^ 

~ 4/15 1/3 1/10 2/15 ■ 1 j 

Pm Pm r m ' 

Here S m is the surface density of planetesimals, M is the protoplanet mass, m is the planetesimal mass, pu and p m 
are the densities of a protoplanet and a planetesimal, Cd is a dimcnsionlcss drag coefficient ~ 1 for km-sized or larger 
planetesimals, and b is the spacing between adjacent protoplanets in units of their Hill radii. An equilibrium between 
mutual gravitational scattering on the one hand and dynamical friction from the planetesimals on the other keeps 
6-10 (Kokubo & Ida 1998). 

The derivation of Equation 2 assumes for simplicity that planetesimal random velocities are always at their equi- 
librium value. Tanaka & Ida (1999) show that rapid migration can significantly boost the planetesimal accretion rate 
of a protoplanet by allowing it to accrete planetesimals with low random velocities existing further inward in the 
disk. However, this assumes that the inner planetesimal disk is pristine. In reality, the shorter accretion timescale at 
smaller radius means that inward migration delivers a protoplanet into a region that has already been gravitationally 
"pre-stirred" by earlier protoplanets of comparable size. Thus we consider our simplification reasonable, and leave 
time-dependent calculation of planetesimal random velocities as a future refinement to the model. More importantly, 
in addition to pre-stirring, these earlier protoplanets have pre-depleted the inner disk. To capture this effect, we track 
the global evolution of the planetesimal disk (Equation 6 below). 

Another assumption which goes into Equation 2 is that protoplanets do not open gaps in the planetesimal disk. 
Since the planetesimal disk has a much lower effective viscosity than the gas disk, a single isolated protoplanet will 
open a gap early on in its growth, unless its migration time is less than its gap-opening time (Ward & Hahn 1995; 
Tanaka & Ida 1999). However, when multiple protoplanets exist in close proximity, planetesimal disk gaps cannot 
form (Tanaka k, Ida 1997; Thommes et al. 2003). 

The net rate of change of the protoplanet mass at a given radius r in the disk (the Eulerian derivative of M) is 
then the protoplanet growth rate (the Lagrangian derivative of M) plus the mass rate of change due to migration (the 
advection term): 



8M(r,t) _ DAI 



dt Dt 



dM 
dr 



(4) 



Chambers (2006) points out that the maintenance of a constant spacing in Hill radius between adjacent protoplanets 
implies that protoplanets must merge, at a rate constituting 50% of the mass accretion rate due to planetesimals. This 
result is obtained for protoplanets that are radially fixed; with migration, protoplanet orbits gradually diverge (as can 
be seen by skipping ahead to Figure 6), so that less merging between neighbors is required in order to keep a constant 
Hill spacing. Nevertheless, in neglecting protoplanet-protoplanet mergers, we are underestimating the protoplanet 
growth rate. We leave a detailed examination of this issue to future work. 

As the protoplanets grow, they deplete the surface density of planetesimals in the nearby part of the disk. If a given 
protoplanet accretes from an annulus of width Ar, then the rate of change of total planetesimal mass in that annulus 
is given by 

d ,„ A ^ x DM 
- (27rrAr£ m ) 



dt v Dt 



(5) 



Using Ar = brn , we obtain from this the surface density rate of change in terms of the protoplanet accretion rate 
(since we neglect planetesimal migration, this equation contains no advection term): 



S£ m (r,i) 1 /M»\ 1/3 DM 



(6) 



dt 3 2 / 3 bwr 2 \ M J Dt 

We can then solve the system of partial differential equations (4) and (6) to obtain M(r,i) and £ m (r, t). 

3. THE DISK MODEL 

In order to model the viscous evolution of the gas disk in a simple way, we make use of the similarity solution 
of Lynden-Bcll & Pringle (1974) We assume an a-parameterization of viscosity, v = ac s H, where c s is the sound 
speed and H ps c s /tl is the gas disk scale height, and assume a to be constant. We adopt a disk temperature profile 
T oc r -1 / 2 (i.e. c s cx r~/ 4 ). The similarity solution is then of the form (Hartmann et al. 1998) 

E ( r i) = Md(0) - e -<r/*>)/r d (7) 

Here Md(0) is the total disk mass at t = 0, i?o is the characteristic disk radius at t = (beyond which surface density 
falls off exponentially), and th = t/(R 2 /3i , o) + 1 is the nondimensional time. For r << RoTd, (3 (Equation 1) « 1. 
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We then scale the surface density of the planetesimal disk to that of the gas disk at t = 0, adding in a factor of 0.24 
decrease inside the water evaporation radius, for which we adopt 2.7 AU as in the model of Hayashi (1981): 

S m (r, 0) = 0.019 x 10 [Fc/H] ./T g (r, 0) (8) 

where [Fe/H] is the metallicity index (0=Solar abundance). We choose [Fe/H]=0.25. The factor / introduces the 
surface density jump at the snow line, smoothed over a radial scale of 1 AU: 



/ = 0.24 +(1 -0.24) 





'r-2.7AU\ 


1" 


- tanh 1 




2 \ 


v 1AU ) 


+ 2 



(9) 



For the sound speed, we use c s = 1.4 x lO^r/AlfT 1 / 4 cms" 1 which yields a disk scale height H « 0.047(r/l AU) 5 / 4 AU, 
again as in the model of Hayashi (1981). In addition, we adopt p m = puj = 1.5gcm -3 , b = 10, and m = 10 -12 M ffi 
(which corresponds to r m ~ 1 km). 

4. COMPARISON OF TIMESCALES 

Before we proceed with solving Equations 4 and 6, we can can gain insight into the problem of concurrent protoplanet 
accretion and migration from a simple analysis of the timescales involved. The accretion and migration timescales of 
a protoplanet are 

- M 

taCCr = DM] AC ( } 

and 

r 

^migr — ~ - ( 1 1 ) 

' migr 

Adopting a disk mass — 0.15 M Q and radius Rq — 50 AU, the model gas disk of §3 above has S 9 (5 AU) ~ 770 gem -2 
at t = 0, and a midplane gas volume density p g ~ T, g /2H, where if (5 AU) « 0.07 AU. We can write the accretion and 
migration timescales as 

t accr ~2.4x 10 5 



E s V 2/5 f H \ 2/5 / r \V2 



(sau) 



M 



10M fl 



770 gem- 2 J \0.07r J \5 AU/ V 0-019 x 770 x 10°- 25 gcm- 

1/3 

yrs, 



S ff Y fH/rYf r Y 2 ( M 
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Thus the accretion timescale of a 10 M ffi body at 5 AU starts out « 6.7x as long as its migration timescale, meaning 
that a protoplanet will fall into the star well before it reaches core size. However, the accretion time has a weaker 

2/5. i 

dependence on gas density than does the migration time; iaccr oc T, g while t m - lgT oc Si" . Therefore if the gas density 
is reduced far enough, accretion can win out over migration. In particular, for the above example, reducing S 9 by a 
factor 6.7 -5 / 3 — 0.04 will make the two timescales equal. 

We can also construct a more quantitative estimate. In the absence of migration, we would estimate the protoplanet 
mass at a time t\ by setting t\ = i accr , then using this together with Equations 2 and 10 to solve for M accr (r, ti). 
The upper bound on accretion is the oligarchic isolation mass, i.e. the mass at which a protoplanet has consumed all 
planetesimals within its feeding zone, Ar = bru : 



MUr) = 2r^ 2 -^. (13) 

With the protoplancts migrating, we can estimate an upper bound on accretion by computing the mass M cross (r) at 
which i accr — t m igr- It then remains to pick a limiting time, after which the disk no longer contains enough gas to 
provide a jovian atmosphere. Using the disk model of §3 with a = 10~ 2 and, as above, Md = 0.15 M Q and Rq = 
50 AU, we examine the time evolution (Figure 1). We compute the times at which the disk has dissipated to the 
point that there is a Jupiter mass of gas left inside 30 AU and 100 AU: t30AU = 2.4 Myrs, iiooAU = 5.5 Myrs. For 

— 1/2 

comparison, the total disk mass at these times, given by Md(0)r d , is 35 Mj up and 24 Mj up , respectively. Also, 
the gas accretion rate onto the star — and thus the approximate mass flux anywhere in the inner disk — is given by 

M d = (3/2)M d (0)(^ /i?o) T £ r 3/2 : and is 7 -° x 10_6 M Jup/yr and 2.1 x 10~ 6 M, Iup /yr, respectively. This suggests t 30A \j 
as well as £iooau constitute reasonable upper limits on how fast a core must form in order to still be able to accrete a 
gas giant atmosphere from the disk. 

Figure 2 shows a plot of Afj SO (r) and M cross (r), together with M accr (r-, £30 au) and M accr (r, tioo au)- In the absence 
of migration, our estimate of the protoplanet mass as a function of radius at time t would be 

M(r,t) = min(Af accr M),M iso (r)). 
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Fig. 1. — Time evolution of an a disk having the sound speed profile of §3 and a = 10~ 2 , initial disk mass Af<j(0) = 0.15 Mq, and radius 
Ro = 50 AU. TOP: The total gas mass (solid), the gas mass inside 100 AU (long- dashed), and the gas mass inside 30 AU (short- dashed). 
BOTTOM: The mass flux in the inner disk in units of Mj up /yr (solid) and M0/yr (dashed). 



With migration, we can estimate the protoplanct mass as a function of (initial) orbital radius as 

M(r,t) = mm{M accr (r,t),M iso (r),M crosB (r)). 

For comparison with the subsequent full computation, in which we use 1 AU as the inner boundary, we take for t m ; gr 
the time to migrate to 1 AU, rather than all the way to the star. To estimate the largest protoplanet mass at time t, 
we find the masses at which M accr (r, t) intersects Mi so (r) and M cross (r), and pick the smaller of the two: 

M max (t) = min (M accr (r, t) D M iso (r), Af accr (r, t) n M cross (r)) (14) 

In this way, we estimate the largest protoplanet mass at £30 au and £iooau to be w 10 M® and 20 M®, respectively. 
Thus, we predict that this set of parameters will indeed produce core-sized bodies, and do so while there is still enough 
gas in the disk to plausibly provide a gas giant envelope. 

5. COMPUTATION OF THE MODEL 

5.1. The simplest case: No accretion, no gas disk evolution 

We now proceed with a full solution of Equations (4) and (6). The equations are solved numerically using Mathe- 
matica 5 and MATLAB 7. We take 1 AU as the inner boundary of our computation domain. For comparison, we begin 
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Fig. 2. — Isolation mass Mi ao (r) (thick dashed curve) and "crossover" mass M croBB (r) (thick solid curve) at which t acC r = tmigr, 
together with M accr (r, t^Q au) ffiii dashed curve) and M acC r (r, 1 100 AU ) ffiin soiid curve) for the gas disk in Figure 1, with all other 
parameters as in §3. At t = t30AUi our estimate of the protoplanet mass as a function of original stelloccntric radius is M(r) = 
min(M accr (r, t30Au), M iso (r), M cross (r)); at t = tiooAUi the estimate is M(r) = min(M acC r (r, iioo Au)> M so (r), M croS s (r)) 




r/AU 



Fig. 3. — Protoplanet mass as a function of stelloccntric radius and time for the case of §5.1: A time-invariant gas disk of mass = 0.15 
M©, radius i?o = 50 AU, metallicity [Fe/H]=0.25, and no Type I migration. The ovcrplotted streamlines show examples of protoplanet 
evolution paths; because there is no migration, all are parallel to the t-axis. 
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Fig. 4. — Protoplanct mass as a function of stcllocentric distance and time with Type I migration, in a time-invariant gas disk of 
mass Mtj = 0.15 Mq, radius Rq = 50 AU (§5.2). The overplotted streamlines show examples of protoplanct evolution paths. The initial 
spacing between streamlines is 0.5 AU, and does not represent the actual radial distance between adjacent protoplanets. TOP: Disk with 
[Fe/H]=0.25; BOTTOM: [Fe/H]=0.5 



by neglecting Type I migration altogether, and we fix the gas disk at its t — value. The result of the calculation is 
shown in Figure 3, for a disk mass of Md = 0.15 M©. A largest protoplanct mass of 10 Mq is reached in w 2 Myrs. 
After 10 Myrs, the largest protoplanet mass is nearly 80 Mg, with this maximum occurring between 5 and 10 AU. 
Thus, in the absence of migration, this protostellar disk readily and rapidly produces bodies of more than sufficient 
mass to furnish jovian cores. 

5.2. Adding migration 

Next, we include Type I migration (Equation 1) and repeat the calculation. Results are shown in the top panel 
of Figure 4. The outcome is now dramatically different. The largest protoplanct mass reached is just over 3 M ffi ; 
this occurs at the inner boundary within the first million years. After 10 Myrs, the largest protoplanet mass is just 
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Fig. 5. — Protoplanet mass as a function of stellocentric distance and time with Type I migration, in an evolving a. disk having a = 10~ 2 , 
Mrf = 0.15 Mq, Rq = 50 AU and [Fe/H]=0.25 (§5.3). The gas disk evolution is thus the same as that shown in Figure 1. The overplotted 
streamlines show examples of protoplanet evolution paths. The calculation is carried forward to tiooAU=5.5 Myrs (see §4), and a line is 
drawn showing t30AU=2.4 Myrs. 



under 2 M ffi . The overplotted protoplanet evolution paths show that Type I migration removes a prodigious amount 
of material: By ~ 5 Myrs, all the protoplanets originating inside 20 AU have fallen through the inner boundary at 1 
AU. We repeat the calculation once more, now increasing the metallicity to [Fe/H]=0.5. This makes the planetesimal 
surface density 46 g/cm 2 at 5 AU, about 17 x the value at that radius in the minimum-mass Solar nebula model. 
(Hayashi 1981). The outcome is plotted in the bottom panel of Figure 4. The mass after 10 Myrs is still less than 
3 Mq, though it briefly exceeds 5 M 9 at the inner edge within the first million years. This illustrates how in a 
time-invariant gas disk, even a very high-metallicity one, Type I migration completely overpowers accretion. 



Dropping the metallicity back down to [Fe/H]=0.25, we now perform the calculation with the gas disk evolving in 
time according to Equation 7. As in §4, we chose a = 10~ 2 , Md = 0.15 Mq, and Ro = 50 AU. The result is shown in 
Figure 5. 

The computation is carried forward to t = <iooau = 5-5 Myrs. By this time, the largest protoplanet mass reached 
is just over 10 M®, around 3 AU. Before this, at £30 au = 2.4 Myrs, the maximum protoplanet mass is « 7 M ffi . The 
masses produced are thus comparable to (though somewhat lower than) the estimates from §4; the model succeeds in 
producing bodies large enough to plausibly serve as giant planet cores while there is still enough gas in the disk to 
furnish a jovian-mass atmosphere. 

Looking at the overplotted protoplanet evolution paths, we see that only the protoplanets originating at < 5 AU are 
lost. It is the short growth timescale at small r which dooms them; they become massive too early, whereas the ones 
slightly further out remain small until the gas disk surface density has dropped to a safer level. 

Figure 6 shows another view of the protoplanets' time evolution, projecting the evolution streamlines into the t — r 
and t — M planes. This gives an even clearer view of the range of outcomes: The promising jovian core candidates, 
which attain a mass close to 10 Mq by tiooAUj originate from narrow radial range, about 5 to 6.5 AU. The ones 
starting further in are lost before they grow massive enough, while the ones starting further out simply don't do much 
accreting at all in the time available. Also visible is the way accretion stalls as each successive protoplanet migrates 
into the region of the disk already picked clean of planetesimals by previous protoplanets. Strictly speaking, the 
mass eliminated at the inner edge is not yet lost since our inner computation boundary is at r = 1 AU rather than 
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Fig. 6. — The protoplanet evolution paths of Figure 5, projected into the r vs. t plane (TOP) and the M vs. t plane (BOTTOM). Over 
the time shown (5.5 Myrs), the innermost eight protoplancts are lost at the inner boundary. The next three are not lost and grow to 8 Mj 
or larger. The remaining protoplancts all remain at masses < 5 Mgj, and undergo relatively little migration. 



r = 0. However, since each of the eliminated protoplanets has already ceased accreting, we are not missing anything 
interesting by removing them a bit prematurely. 

6. DISK PROPERTIES AND CORE FORMATION 

We now investigate how the model results change with disk properties. We already know that the disk must dissipate 
in order for a window of opportunity for core formation to open; now we wish to quantify the dependence of the model 
on disk viscosity (as parameterized by a). We also want to examine how the outcome changes with initial disk mass 
and metallicity. We can already anticipate the qualitative form of the latter dependency: Since increasing metallicity 
increases the planetesimal disk mass, thus tipping the balance away from migration and toward accretion, we expect 
that for a given set of disk parameters, increasing [Fc/H] will always produce larger protoplanets within a given 
time. However, dependence on initial gas disk mass is not a priori obvious. Since we take S m cx S s (0) oc Md(0), 

the growth rate at t = goes as dM/dt oc Sg^ 5 S m oc M d (0) 7 / 5 . Meanwhile, the migration rate at t = goes 
as dr/dt oc S g oc M^(0). Thus, increasing M<j(0) boosts the initial advantage of accretion over migration, but as 
illustrated in Figure 4, an initial spike in protoplanet mass can be rapidly erased by migration. Figure 7 repeats the 
computation shown in Figure 2 with four different disk masses, keeping a = 10~ 2 . We see that initially, M max is set 
by the isolation mass Mj so rather than by Af cross , and increases with M^(0). M CTOSS decreases with disk mass while 
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Fig. 7. — The calculation shown in Figure 2, repeated for disk masses of 0.03, 0.05, 0.1 and 0.4 Mq. In each case, a filled circle shows 
the estimate of M max at t = tiooAU ■ 



Mjso increases, until M max is set by M cmss . For even larger disk masses, M max then remains approximately constant 
but occurs at ever-larger r. 

We can rapidly scan a large parameter space using the estimation scheme of §4. For each set of parameters, we 
can calculate M max (Equation 14). We proceed by varying Md(0) and a, holding all other parameters fixed at the 
values used for §5.3. We consider disk masses from 0.01 to 0.4 Mq, and a from 10 -4 to 10 _1 . We repeat this for 
[Fc/H]=-0.25, 0, 0.25 and 0.5. The results are shown in Figure 8. As expected, M max increases with metallicity. Also, 
we see that for a given value of the disk mass, there is indeed an "optimal" value of a, = abcst, which produces the 
largest M max . If a <C abest, the disk does not dissipate fast enough and significant migration takes place, resulting in 
an outcome similar to those depicted in Figure 4. On the other hand, if a 3> o^cst the gas disk dissipates too fast, so 
that by t = tiooAUi little accretion has occurred. The value of cibest increases with disk mass and metallicity. Finally, 
we see the dependence on disk mass indicated by Figure 7: for a given value of a, M max initially increases with Md(0), 
then levels off. The lower a, the lower the disk mass at which M max stops increasing. 

We now calculate a grid of full models as in §5.3 for the same parameter range as above, in each case extracting 
Af max at t = iiooAu; the results are shown in Fig 9. Relative to Figure 8, M max is less strongly peaked toward 
large disk mass and high viscosity. At low viscosity, rather than leveling off, M max increases monotonically with 
initial disk mass for a given a. Our simple estimate nevertheless fares relatively well in qualitatively reproducing the 
Mjax — Mrf(0) — a — [Fe/H] relationship, and over the plotted range gets the maxima of the M max surfaces right within 
a factor of < 1.5. 

From Figure 9, we see that a sub-Solar metallicity disk, [Fe/H]=-0.25, falls just short of producing M max = 10 Mq, 
even with an initial disk mass of Md(0) = 0.4 M©. A Solar-mctallicity disk produces M max = 10 M© if it starts out 
with Md(0) > 0.3 M©, provided also a ~ 10~ 2 . With a super-Solar metallicity of [Fe/H]=0.25, a disk mass slightly 
above 0.1 M© can produce M max = 10 M© if, again, a ~ 10~ 2 . Finally, with a very high metallicity of [Fe/H]=0.5, a 
disk mass < 0.1 M© suffices to make M max = 10 M©, if also 10~ 3 < a < 10~ 2 . 

7. DISCUSSION 

We have demonstrated that bodies ~ 10 M© in mass, large enough to initiate the formation of gas giants by nucleated 
instability, can grow in a viscously evolving gas disk even if they are subject to Type I migration. Though the accretion 
timescale is significantly longer than the migration timcscale in a gas disk that is young and massive, Md ~ 10 _1 M©, 
the times become comparable when Md ~ 10~ 2 — 10 -3 M©; in other words, when there are only a few Jupiter masses 
of gas left in the disk. At this point, a protoplanet has a chance of growing to core mass before spiraling into the star. 
In this picture, then, the formation of giant planet cores does not take place until a window of opportunity opens in 
the latter part of the gas disk's lifetime. 

At the same time, though, our findings suggest that an era of successful core formation is by no means inevitable in 
a protoplanetary disk. The accretion of 10 M© or larger bodies necessitates an initial disk mass Md(0) > 0.1 M© at 



10 



Thommes & Murray 





Fig. 8. — The largest protoplanet mass, M max (M®), at the time the gas disk mass is 1 Mj up inside 100 AU, using the estimate of §4, 
for [Fe/H]=-0.25, 0, 0.25 and 0.5. 



Solar metallicity, and not much less than 0.1 Mg at [Fe/H]=0.5. A lower-metallicity disk having [Fe/H]=-0.25 does 
not produce 10 M® bodies unless it is very massive, Md(0) > 0.4 Mq. Furthermore, at a given disk mass the largest 
protoplanets occur only at a particular viscosity. Too much lower, and the disk "overstays its welcome" , removing too 
much mass in abortive cores from the disk. Too much higher, and the disk dissipates too rapidly, without core accretion 
having had a chance to get very far. However, assuming disks to be adequately described by an alpha parameterization 
of viscosity, the values we obtain for «bost, ~ 10 -2 and as low as 10 -3 , agree well with values of a obtained from fits 
to observed disk lifetimes and accretion rates (Hartmann et al. 1998). We have not considered disk photoevaporation 
(Shu et al. 1993; Matsuyama et al. 2003). In a disk where this exists as an additional sink of disk gas, abest will be 
correspondingly lower than it would be for pure accretion onto the star. We defer a detailed examination to future 
work. 

It has been suggested (e.g. Trilling et al. 2002) that protostcllar disks produce numerous gas giants, but that inward 
migration in an annular gap (Type II migration) removes many of them as the disk accretes onto the star, leaving 
only the last few survivors to populate the mature system. The model described here is analogous to this picture, 
except that the culling of giant planets occurs already at the core stage. In fact, our analysis suggests that full-sized 
gap-opening gas giants do not have a chance to grow until the gas disk is quite tenuous. Since the runaway accretion 
of the gas envelope can proceed at a rate > 10~ 5 Mj up /yr (Pollack et al. 1996), while by this time the disk accretion 
rate is ~ 10 -6 — 10~ 5 Mj up /yr (§4), a significant fraction of the disk's accretion likely goes into feeding the planet 
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[Fe/H]=-0.25 [Fe/H]=0.0 





Fig. 9. — The largest protoplanet mass, M max , at the time the gas disk mass is 1 Mj up inside 100 AU, computed using the full model, 
for [Fe/H]=-0.25,0, 0.25, 0.5. 



rather than pushing it inward. Once the planet grows massive enough (~1 Mj up ) to open a deep gap and enter the 
Type II regime, its mass is comparable to that of the remaining disk, so that significant Type II migration is unlikely 
to occur. Therefore in our scenario it becomes Type I rather than Type II migration which plays the primary role in 
determining the radial distribution of gas giants. 

The above argument notwithstanding, a detailed examination of the endgame of giant planet formation in a tenuous 
gas disk is an important area for future work. Is our adopted criterion of a Jupiter mass inside 100 AU a sufficient 
condition for gas giant formation? Pollack et al. (1996) find that even after a core grows to ~ 10 M ffi , a plateau of slow 
gas accretion can exist for many millions of years before runaway gas accretion supplies the massive jovian atmosphere. 
If this is the case, lots of extra gas will be required in the disk. However, these model used the interstellar value for 
grain opacity; reducing this to account for settling and coagulation of dust in a protostcllar disk can substantially 
reduce the length of the plateau, as shown by Hubickyj et al. (2005). The authors also demonstrate that cutting off 
planctcsimal accretion onto a core speeds the onset of runaway gas accretion, by reducing the rate of energy deposition 
into the envelope. In fact, migration does exactly that, as illustrated in Figure 6: As protoplanets move into parts of 
the disk already picked clean of planetesimals by their inner neighbors, they cease accreting. 

Once the first gap-opening planet does form, it will act as a barrier in the disk to faster migrators originating from 
outside its orbit (Thommes 2005), thus the formation of subsequent planets — including gas giants if there is enough 
gas left in the disk — becomes easier. Before that, however, a significant total mass in failed cores may migrate to 
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the central star. Recent discoveries of short-period exoplanets with a minimum mass comparable to that of Neptune 
orbiting GJ 436 (Butler ct al. 2004) and 55 Cancri (McArthur et al. 2004) may represent such would-bc cores, which 
survived either through luck (the gas disk was removed just in time) or through some migration-stopping mechanism 
near the star. 

The observed correlation between host star metallicity and planet occurrence rate (Gonzalez 1997; Fischer & Valenti 
2005) is seen as evidence that gas giants begin their existence as solid cores, since more solid material allows a disk 
to form larger solid bodies, and to do so faster. The model of concurrent core migration and accretion developed here 
provides an even stronger theoretical basis for a planet-metallicity correlation: The higher the ratio of solids in the 
disk, the better accretion does relative to migration, all other things being equal. 

Finally, as alluded to in §1, we have considered here what is essentially the worst-case scenario of Type I migration. 
Details of the disk structure, not considered in our simple model, may act to significantly reduce migration. The 
MHD turbulence providing the disk's viscosity may simultaneously produce enough fluctations in the gas density (and 
thus in the disk torques) to give Type I migration the character of a random walk in radius (Laughlin et al. 2004; 
Nelson & Papaloizou 2004), thus enabling at least a tail of the protoplanet population to survive even in a massive 
disk. Persistent rather than random jumps in disk properties may do even better. Mcnou & Goodman (2004) show 
that at the locations of opacity transitions, Type I migration may slow down by an order of magnitude. The density 
jump associated with the outer edge of a dead zone (Gammie 1996) may stop Type I migration altogether (Matsumura, 
Thommes & Pudritz, in preparation). Help may also come on the other side of the migration-accretion contest. As 
mentioned in §2, our simple accretion rate estimate neglects (i) planetesimal migration by gas drag, (ii) planetesimal 
fragmentation, (iii) the growth of gas envelopes on protoplanets, (iv) enhanced accretion before planetesimal random 
velocities catch up to their cquilbrium values, and (v) the contribution to protoplanet growth by mergers among 
protoplanets. Effects (ii)-(v) enhance accretion, and although (i) by itself reduces accretion efficiency (Thommes et al. 
2003), this is largely offset (Inaba et al. 2003; Chambers 2006) by (ii) and (iii): smaller planetesimals experience 
stronger gas drag, which increases their migration, but also increases the efficiency with which they are captured 
by protoplanet atmospheres. Thus, our estimate of the growth rate is likely conservative. Also, we have assumed 
that the ratio between the solid and gaseous components throughout the protostellar disk at t = is just given 
by the metallicity; in reality, large local density enhancements in the solids may develop (Stevenson & Lunine 1988; 
Cuzzi & Zahnle 2004), giving a correspondingly large boost to core accretion at those locations. 

In short, there is an abundance of candidates for countering the effect of Type I migration. Our purpose here 
has been to demonstrate that even in the absence of any such mechanisms, Type I migration is far from being an 
insurmountable obstacle to the formation of gas giants by core accretion. 

8. SUMMARY 

The formation of giant planets by core accretion requires the growth of solids bodies ~ 10 Mg in mass. Bodies of 
this mass are well below the threshold for gap formation, and will thus migrate inward relative to the gas due to the 
imbalance between inner and outer planet-disk torques (Type I migration). In a gas disk with mass of a few x 10 -2 M Q 
or greater, the Type I migration timescale of core-mass bodies is shorter than their formation time, so that growing 
protoplanets plunge into the star before they have a chance to become cores. We have demonstrated here that cores 
can form nevertheless, and that it is the "late bloomers" which are actually favored. Protoplanets which grow large 
early suffer strong migration and are lost. Those which grow later, however, do so in a more tenuous gas disk, with a 
correspondingly longer migration time. A window of opportunity can thus open between the time when growth and 
migration timescales become comparable, and the time when the disk no longer contains enough gas to furnish a jovian 
atmosphere. The size of this window increases with metallicity and initial disk mass. 

We thank Doug Lin for numerous stimulating and informative discussions on this topic, Sarah Nickerson for assistance 
with the MATLAB version of the computations, and the referee for valuable comments and suggestions. This work 
was supported by NSERC of Canada. 
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